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Abstract 


The A-contractivity of Runge-Kutta methods with respect to an inner- 
product norm, was investigated thoroughly by Butcher and Burrage 
(who used the term B-stability) . Their theory is here extended to 
contractivity in a region bounded by a circle through the origin. 

The largest possible circle is calculated for many known explicit 
Runge-Kutta methods. As a rule it is considerab7y smallex* than the 
stability region, and in several cases it degenerates to a point. 

An explicit Runge-Kutta method cannot be contractive in any circle of 
this class, if it is more than fourth order accurate. The practical 
relevance of this analysis is not yet quite clear. 



1. Introduction 


We investigate the contractivity of Runge-Kutta methods when applied to 
nonlinear differential equations. While stability of a method is con- 
cerned with the bovmdedness of the numerical result, contractivity 
requests that the difference of any two numerical solutions, computed 
with the same stepsize, does not grow in a certain norm. For one-step 
methods and the natural norm, given by the differential equation, both 
concepts are identical if the differential equation is linear with 
constant coefficients. In the other cases contractivity is a stronger 
requirement. 

For linear multistep methods contractivity has been introduced by 
Dahlquist [UJ, where it was called G-stability. G stands for a positive 
definite matrix which is method dependent and is used to define a norm 
in the space of numerical solutions. Nevanlinna and Liniger [10] have 
treated contractivity of linear multistep methods using method 
independent noms, such as the maximum norm. Butcher [3] introduced 
B-stability which is contractivity for nonlinear, autonomous contractive 
differential equations using the natural norm. In (I] similar contractivity 
concepts have been discussed, namely AN-stability for nonautonomous linear 
and BN-stability for nonaittonomuous nonlinecu* systems. These concepts 
reduce to A-stability in the linear constant coefficient case and are 
thus only reasonable for implicit methods. We extend the contractivity 
concept for Runge-Kutta methods in such a way that explicit methods are 
included too. We will be using the natural norm in contrast to [2] where an 
idea similar to Dahlquist? s G-stability is introduced. In all these concepts 
one requests a certain monotonicity condition for the differential equation 
In the present article this condition is given in (2.9) • Then it is shown 
that the numerical method when applied to such a differential equation 
is contractive for either arbitrary or special choices of the stepsize h. 

In the remaining part of this section we give an outline of the article. 

In Section 2 basic notations and definitions are given. In particuleu* 
the monotonicity condition for the nonlinear differential, equations and 
the concept of contractivity sure described. In Section 3 the r-circle 
contractivity is introduced. If a method is r-circle contractive then 
the stability region contains the interior or exterior of a disk of 
radius |r| which is tangential to the imaginary axis at the origin. 

However, the converse is not true, i.e. there are methods whose stability 
region contains a disk of radius r with the origin on the boundary 
which are not r-cirele contractive. We then give purely algebraic necessary 
and sufficient condition in terms of the coefficients for a method to be 



- 2 - 


r-circle contractive. An alcoritba is given vhich enables one to coaqpute 
r for any given explicit or is^l:' sit Runge-Kutta method. It is natural to 
introduce the concept of reducible methods. 

An m-*stage Bunge-Kutta method is reducible if there exists an m*-stage 
Ihinge-Kutta method vith m* < m and both methods give identical results on 
any computer >diich carries out additions of 0 and multiplications by 0 
without round-off errors. It is then shown that for irreducible r-circle 
contractive methods 1/r is a continuous function of the coefficients of 
the method and that this is not the case if one admits reducible methods. 
Further confluent methods are introduced. A method is called confluent if 
at least two of the row sums of the coefficient matrix A are equal. It is 
then shown that to any confluent method* which is r-circle contractive and 
to any e > 0 there exists a nonconfluent method which is r* circle contractive 
and 1 1/r - 1/r* I < e. In Section U we show that one has numerical contractivity 
for nonlinear differential equations if the swthod is r-circle contractive* 
if the differential eqtiation satisfies the monotonicity condition (2.9) and 
if h is chosen appropriately.* In Section 5 ve show that for ui explicit 
r-circle contractive method one has r<m* where m is the number of stages. 
This result is sharp. Further if r is negative then the p is 1 error order 
and r < 1/2c where c is the error constant of the method. Finally* we list 
the value of r for many of the well known explicit Runge-Kutta methods. 


* It is shown that an explicit circle contractive sMthod cannot have an 
error order exceeding 4. 
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2. The methods and the test equation 


For solving initial value problons 

(2.1) y'(t) ■ f(t,y(t)), y(0) given, y,f € R* or C® , 

ve consider m stage Runge'-Kutta methods. Let h>0 be the step site, t^*nh 
and y^ is the numerical approximation to the exact solution y(t^). The 
numerical solution y^^^ at t^_^^ ■ t^ + h is canputed as 


( 2 . 2 ) 


vhere 




(2.3) 


Y. * y ^ ^ ^ l*l92)aea)lBa 

^ n u j j 


We always request 


(2.Ua) 


m 

I ■ 1 


(2.»*b) 


m 

y a. . . 

j-1 


Observe that by (2.Ua) the method heus an error order of at least one. 

(2.Ub) is not necessary for a method to be convergent, see [ll]. 

However, it is convenient in notation to have (2.Ub) and practicedly all 

known methods satisfy (2.Ub). Moreover, the extension of the present 

results to methods without (2.Ub) is trivied. If the matrix A > (a. .} is 

^ J 

strictly lower triangvilar then the method is ccdled explicit otherwise 
implicit. We cedi a method nemoonfluent if all c. are distinct and 
confluent otherwise. For compactness in notation we introduce the vectors 
Y,F^(Y)6IR“® or c“® and defined by 


(2.5) 


/fCV',".!,) 

VY)!-/ f(t„«jh.Y5) 


i n m m 4 


We shall simplify the notation by the use of the K' necker*^roduet symbol V, 
see [ 5, p.1l6]. In order to avoid parentheses ve assume that • has higher 

priority than ordinary matrix multiplication. Let I be the sxs identity 

. T / * 

L.atrix, and let b ■ (b^,bg,...,b ). Ihen (2.2) reads 


( 2 . 6 ) 




and (2.3) takes the form 


(2.7) 


Y ■ l*y„ + hA»I F (Y) . 
n s a 
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The aim of this article is to investigate under vhat conditions any 
tvo numericed solutions ly > * f which are ccmputed 

with the same h vill satisfy the inequalities* 

(2.8) ’*yn+l“*n+1** - *’^n“*n*’* “"O*^***** 

We assume here that Hull :> ^u^u^^ where is an inner product defined 

on or C^. Note that in contrast to G-stahility [U] and the nonlinear 
stability in [2] the norm does not depend on the method used but only on 
the differential equation treated. We tadk of numrioal oontraetivity if 

(2.8) is satisfied. The main purpose of this article is to show numericed 
contractivity. To do this we need to impose conditions on the differential 
equations and on the methods. The condition on the method is the r-eircle 
contractivity which is treated in Section 3. For the differential equation 
we request the monotonicity condition 

(2.9) Re^f(t,y) - f(t,z), y-z) < -a llf(t,y)- f(t,z)l? vy,z€R®orC®. 

In Section h ve shall shov that if a*r and the stepsize h satisfy the 
inequality (U.2) then one has numerical contractivity. To clarify the 
condition (2.9) we observe that for a linear equation y* ■ Xy condition 

(2.9) becomes 


( 2 . 10 ) 


Red <t-aX)/X < 0 . 


Thus if ve introduce the generalized disks 


(2.11) D(r) 


(X€C| 

|X + r| <r) 

if r> 0 

{X€5| 

ReX a 0} 

if r««» 

(X€C| 

|X+r| k -r) 

if r< 0 


then (2.10) is equivalent to X€D(l/2o). If a>0 then (2.9) implies that 
for tvo solutions y(t) and z(t) of (2.1) one has 


(2.12) ^lly(t)~z(t)lP<0 for all t. 


Further observe that a is not invariant against sealing. Let y(t) be a 
solution of (2.1) and define z(t):« y(rt). Then z(t) is a solution of 
the scaled system 

*'(t) ■ g(t,z) 


where 


g(t*z):* Tf(tT,y) 



If (2.1) satisfies (2.9) with a“o^ then g satisfies (2.9) with o«a “to.. 

* g I 

Moreover (2.9) with a>0 implies that f is Lipschitz continuous with 1/a 
as Lipschitz constant, for one has 

llf(t,y) - f(t,z)ll lly-zll > Re^- f(t,y) + f(t,z),y-z) 

> aRe{- f(t,y) + f(t,z), - f(t,y) + f(t,z)^ 

-Re^f(t,y)- f(t,z),y-z + of(t,y) - af(t,z)) >ollf(t,y) - f(t,z)l|2 . 


Here we have used Schwarz's inequcdity and (2.9). 
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3. The r-circle contractivlty 

In this section ve define r-circle contract ivity. In order to motivate 
this definition we consider the scalar test equation 


(3.1) y* ■ X(t)y(t), X(t)€C. 

If one applies (2.6)* (2.7) to (3.1) the numbers 

(3.2) ■ hX(t^+hCj^), i«1,2,...,m 

&nd c B (c^ tCgt. . • tC^) are needed. Assume that (3.1) satisfies the 


monotcnicity condition (2.9) then Cj^6D(r) with r“h/2o. If the c. are 
distinct then one can choose any m complex numbers €D(r) and find a 
smooth X(t) such that (3.2) holds. Applying (2.6), (2.7) to (3.1) leads 


to 


(3.3) 


where 


(3.M 

rt(t) ■ 1. Ad -Az)->1 

Bl 

with 


(3.5) 

Z ■ diag( g • • • t 


see [l]. Clearly we have numerical contractivity if |k(c)| < 1. This leads 
to the 


Definition 1 . a Runge-Kutta method is called T-oiTol« aontraetive if D(r) 
is the largest generedized disk with reO and 

(3.6) Ik(c)| < 1 for all c€D®(r). 


A method is called oirol9 oontraoHvt if (3.6) holds for some reO. 


Note that for a confluent method applied to (3.1) one never has (.#<;• 

^ J 

if c. *c.. Nevertheless we request (3.6). One reason for this is, as we 

* J 

shall see at the end of this section, that with the present definition 

1/r is a continuous function of the coefficients a., and b. if the method 

J 

is irreducible. Clearly D(r)cs, where S is the stability region of the 
method, given by 


8 - {y€c| |K(vl)| S 1} . 

Following Burrage and Butcher [l] ve introduce the axatrix. 


(3.7) 

Q ■ BA + A^B-bb* • (iLiY “ 

\^JA»1,j-1 

where 

(3.8) 

B ■ « 
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Theorem 3.1 . 

(3.9) 

and p • - 1 /r 

( 3 . 10 ) 


A Runge-Kutta method is r-circle contractive if and only if 
b. 2 0 for j ■ 1 ,2,... ,m 

V 

is the largest number such that 

i p Bw for all w € . 


Proof . According to Corollary U.3, the conditions (3.9) and (3.10) with an 
arbitrary p* imply that (3.6) holds with r* ■ -1/p' if p*#0 and r'«« if 
p' *0. We then only need the converse result, namely 


Lenina 3.2 . Assume (3.6) holds for some r'eo, r* may be infinite. Then 
(3.9) and (3.10) hold for p'«-1/r* if r' is finite and p'«0 otherwise. 

To show Lemma 3.2 we need the following lemma of Burrage and Butcher [l]. 


Lenina 3.3. Let Z be such that I_-AZ in nonsingular and let 

■ m 


(3.11) 
Then 

( 3 . 12 ) 


u ■ (I - AZ)-ll. 


|K( 


a ® ^ . 

c)P - 1 - 2 [ b.|u.|2Re?^ - I Ci Ui CiU, . 

i«1 i,j«1 ^ ^ J 


Proof of Lenina 3.2 . Assume that for some r' one has (3.6), that is 
(3.13) |K(c)| 2 - 1 s 0 for all c6D®(r'). 

To prove b. > 0, assume on the contrary that b. < 0 for some i. Choose ■ 0 

J ^ BI ^ 

for jai and C£*-e. For c>0 sufficiently small one has c€Ir(r'). By 
(3.11), 

(3.11+) Uj ■ 1 + ^j(e), where |♦j(c)|■*0 as c-*0,j ■ 1 ,2,... ,m . 

The right hand side of (3.12) becomes 
(3.15) -2b^e + ck(c) 

with |k(e)|-*0 as e-*0. (3.15) ii positive for e sufficiently small. This 
contradicts (3.13). 
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T/ 1 

In order to show that v (Q+ 7 rB)v is nonnegative we assume the contrary. 

T mM " 

Let w ■ (w, ,w.,. . . ,v_) €« he such that 
1 C n 

(3.15) I < 0 . 

i,j«1 ^ J i»1 ^ ^ 

i(p.e 


Let (p. ■ w./r' and ?. ■ - r' 


+ r*e ^ ■ i w. 


'j y~ — ’j - - - - "j ® ®y construction 

C ■ (i»1**) holds 

again. We substitute (j in the right hand side of (3.12) and find 

(3.16) |K(c)| 2- 1 • b. W.2- J qii w. w.)c2 + e2k.(c) 

' i«1 ^ i,j«1 ^ ^ 


with |k^(e)|-*0 as c-*0. Hence (3.l6) gives a contradiction to (3.13) for e 
sufficiently small. Thus (3.10) holds for p'>~1/r'. This proves Lemma 3.2 
and Theorem 3.1. • 


Remark . From Theorem 3.1 follows easily that an algebraically stable method 
in the sense of Burrage and Butcher [l] is r-circle contractive with a non- 
positive r. 

In order to describe the situation where sane of the b. are equal to zero 

J 

it is convenient to introduce the 

Definition Z , An m-stage Runge-Kutta method is called reducible if there 
exist two set'.s S and T such that SflT*^, SUTB{1,2,...,m) and 

(3.17) \ ° • 

(3.18) a^j^ ■ 0 if j€T and k€S. 

The method is called irreducible if it is not reducible. 

This definition says that the stages with index in S don't have an influence 
on the final outcome of the integration provided m\iltiplioations by 0 
and additions of 0 are performed exactly. If the method is reducible it is 
equivalent to the m* -stage Runge-Kutta Method which consists of the stages 
with index in T only. Hence m* is the nxasber of eleawnts in T and m* < m. 

We study now Theorem 3.1 for r-cirele contractive methods with some 
Let 8 and T be such that BUT »{1,2,...,m) and 

(3.19) 

( 3 . 20 ) ^ 

By (3.7) • ■ 0 for k€8. Hence for Q-pB to be nonnegative definite 
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it is necessary that 
(3.21) q, 


■kj 


for all k^S. 


Since <ij^j ■ when one finds that (3.21) is satisfied if and 

only if 


( 3 . 22 ) 


a^j^ ■ 0 whenever j € T and k € S . 


Thus (3.19)« (3.20) and (3 >22) imply that the method is reducible. We 
have therefore shown the 


Corollary 3.4 . An irreducible Runge-Kutta method is r-circle contractive 
if and only if 

(3.23) ^ ® 

and p ■ -1/r is the largest number such that 

(3.2U) w’^Qw 2 pw^'bw for all w€IR“. 


Let t/} be the set of all irreducible circle contractive Runge-Kutta methods . 
Hence, by Corollary 3.^, a Runge-Kutta method is ine/^if and only if all 
b. are positive. The methods in are the ones which interest us. If a 
method is not inc/^ic is either not circle contractive or it is reducible 
and after deleting the irrelevant stages one has a member of c^. 

In the following we shall coomute rfor a given method in«^ Since all b. 

U V« U (4 ^ 

are positive it follows that ■ diag(b, , b?,..., b “) is nonsingular. 

Using the transformation B^w «v reduces (3.2U) to 
(3.25) v^B’^^^og'^y j pyTy fQj. ^2^2 v6l^. 


Let v^,V 2 i...,v^ be the eigenvalues of the real and symmetric matrix B'^QB**^ . 

Hence the largest p for which (3*25) holds is p_. ■ min v. and thus 

i-1,2 t " 

by Corollary 3.^ one has 


( 3 . 26 ) 


f 

• if sdn V. ■ 0 

i*1,2,...,m ‘ 

i^1 ,2 , . « . ,m 


Clearly the set««^ is open and p .. ■ min v. ■ 1/r is m continuous 

““ i-1,...,m^ 

function of the coefficients a^j and b^ of the methods. However, if some 
of the bj tend to zero the following possibilities can occur. Either the 
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liaiting method is no longer r-cirele eonti motive, see for example Heun's 
method ($.10), or else it must become reducible. In the latter case r may 
or may not depend continuously on bj, as the following example shova. 


Example 1 . Let 


A ■ 


''0 

^0 


o'^ 

oJ 


b^ - (l-e,c) . 

Clearly (2.U) and (3.9) are satiafied for all c€[0,l]. For e€(0,1) 
ve find 


B 




/ c-1 




The eigenvalues are 

'*1 2 * 2 ^ * ^(2a+1 )2 - 8oc) . 


Hence 


limp . (e) 
c-K)+ 


-1 if o>-'^ 

2o if a<-'j . 


However, if e • 0 then the method is reducible and cui be reduced to Euler's 
method with 

A ■ (0) 
b^ - ( 1 ) 
and 

Hence one has a discontinuity on iff if a<-'^, and if ai-'^ ^min^^^ 
continuou in (0,1 ). 


Hote that similar discontinu ties of can occur as some bj tend to 
zero even if one restricts oneself to the class of explicit methods. 

Observe that the set C of confluent stethode int/i is a surface in of 
lower disMnsion. Thus by continuity of 1/r as a function of a^^ and bj 
any confluent r-circle contractive method in«^ cat. be approximated by a 
non confluent r* -circle contractive method such that 1/r is as close 
to 1/r* as one wishes. This property would not hold if ve had 
replaced (3*6) by 


(3.27) 


|k(c)| si for all c€nP(r)nV 
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where 


V • {5€C?“| c. ■ c. whenever e. ■ c,} 

J 1 J 


as we can see in the following 


Example 2 . 
given by 


Consider the classical 3-stage Nystron method of order 3 


A ■ 



0 

0 

2 

3 






I 





see [9] p.U8. If one computes r using the above algorithm one obtains 
r«0. 92668857. If we haa used (3.2?) instead of (3.6) in the 
definition of r-eircle contractivity one would have found **g®3. However, 
for ~ e sufficiently smeU.1 (3.6) and (3.27) are identical. Thus using 

(3.27) instead of (3.6) would have resulted in an r which does not dei>end 
continuously on the coefficients of the method. This is one reason for 
choosing (3.6) rather than (3.27). The main reason, however, is the 
Theorem U.1 of the next section. 

Condition (3.27) was used by Burrage and Butcher [I]. A similar condition 
(with intervals on the negative real axis) was used as early as in 1957 
by Liniger [l3l. 
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4. Nonlinear contractivlty 

Theorem 4.1 . Assume the differential equation satisfies the monotonicity 
condition (2.9) and the Runge-Kutta method is r-circle contractive. Then 
two numerical solutions and computed using the same stepsize h > 0 
satisfy 

*’^n+1 “ ®n+1** “ **^n" *n** “*0»’»2.*** 

provided 

h/r<2o if ra« 

( 1 » 2 ) 

a > 0 and h arbitrary if r = «» . 


Proof . First we observe that it is enough to show (U.l) for n»0 only. 

Subtracting from (2.6) the corresponding equation for the solution 

\ gives 

i ‘^■'11*0,1 ,. . . 


X = X- + hb • I P 
1 0 s 


(U.3) 

where we have used the abbreviations 


xo®yo"2o» F * Fq^y) ~Po(2) 

and Z € If® or C“® ia given by 



In a similar fashion one obtains from (2.7) the equation, 

(**.5) X»l*Xrt + hA«IF, 

where X ^ y - Z. It is enough to show that 

(U.6) llxil|2- IIxqIP s 0. 

By (U.3), 

(U.7) lIx^lP - IIxqIP ■ h2Re(xQ,b’^ • I^f) + h^ Hb’^ • I^FlP . 

The first term on the right hand side can be simplified if we introduce 
the following product t , 1 in if* or C?“. Let U,V € if® or <f® be given tgr 
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(U.9) (xq. « [1«Xq.F1 . 

In order to show (U.6) we need an upper hound for Retl«XQ,P]. The following 
lemma is an easy consequence of (2.9) and the definition (U.8). 


Lemma 4.2. Assvmie b.>0 for j « 1,2,... «m and that the monotonicity 
■ J 

condition (2.9) holds. Then 

(U.10) Re[P,X + oF] < 0 . 

Eliminating X frcsn (U.10) using (U.5) leads to 
(U.11) Re[P,l®XQ] s -hRe(P,A»IgP + ^P] . 

Using (U.9) and (U.11) in (U.7) gives 
(U.12) lIxilP- !!xolP < -h2ReP(P) 

where 

(U.13) P(P) - 2l(A®I )P,Pl- Ilh’^wlPlP. 

s n ns s 

Observe that P(P) is a quadratic form in F and it remains to show that its 
real part is nonnegative. Let G € 1^® or d?®® be written as 



where € R® or C®. Hence 



-1U. 


itep(o) • X 






Thus by (3.10) ReP(O) is nonnegative if -2o/h < p ■ -1/r if r*«. 

If then a has to be nonnegative and h is arbitrary. This completes 
the proof of Theorem U.1. 


Note that (U.2) also shovs that the scheme is numerically contractive in some 
cases vhen the differential system is not so. namely if a<0 and r< 0. 


Corollary 4.3 . Assume that 

(U.lU) b. 2: 0 for j«1,2,...,m, 

J 

and 

(1».15) w^Qw i p'w^Bw for all w€l^ . 

Then one has 


(U.16) 


|K(c)| S 1 for all c€D"(r») 


where 


r' 


• if p' ■ 0 
-1/p' otherwise. 


Proof. Let c ■ € D^(r'), Z « diag(c,.... Let h» 1, s« 1 

and P«ZX where X € Iben {k,3) and (U.5) become 


(•♦.IT) 

xi ■ xo ♦ b^ZX 

and 


(J*.18) 

X ■ XqI -I- AZX . 

Thus 


(U.19) 

Xi ■ K(c)xo 


and we have proved the corollary provided (U.6) holds. This is. however, 
shown in exactly the same way as in the proof of Theorem U.l. Just observe 
that c€l^“(r') implies RetZX.X-f ezx] s 0 for a»l/2r* and that (1*.15) with 
p* > -1/r* implies ReP(0)i0. 
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Theorem 4.4 . An irreducible Runge-Kutta method that is more than fourth 
order accurate » cannot be circle contractive (with respect to the norms 
considered in this paper), unless 

1 “ 

(U.20) ■ I i*1,2,....m. 

1 ij J 

These conditions cannot be satisfied, if the method is explicit. 

Proof . We first prove that (U.20) cannot hold for explicit methods. Since 
for an explicit method c. ■ 0 and a. . ■ 0 for j > i, it would follow from (U.20) 

I X j 

by induction that c^ * 0 for i ■ 1 ,2,3,. . . ,m, which is possible for first order 
accurate methods only. 

The other statements follow from the positivity of the b. (Corollary 3.1») 

J 

and a lemmas that John C Butcher pointed out to us. 


Lertina 4.5 . 
accurate • 

(U.21) 


(Butcher [lU]). If a R\uige-Kutta method is more than fourth order 
then 


f b. (? a..c.-ic.2V = 0. 


Proof. The left hand side of (U.21) can be written. 



These sums are associated with the rooted trees of order 5 named, respectively, 
ti 2 , tg and tg in Butcher's euLgebraic theory, and, by tl5» Table 9.3] their 
vcdues are, for any method that is more than fourth order accurate, respect- 
ively, 1/,’0, 1/10 and 1/5. Hence 


L 




1 1 


TT’5 


0 . 


This proves the lemma and the theorem. 


The conditions (U.20) cmd more general conditions of a similcu: type are often 
x'equired for implicit Runge-Kutta methods, see e.g. [3l. 

See C.XSO the comments at the end of the paper. 
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5. Methods with optimal r and examples 

Given an r-circle contractive Runge-Kutta Method, let D(r^) he the largeot 

generalized disk of the form (2.11) in ita ctahiliiy region H. Thon one has 

D(r)cD(r^). The following two excunples show that D(r) may be a proper 

subset of D(r ). 

s 


Example 3 . The e-method is given by 



b'^ ■ (0 1-e) 


or 


'n+1 






» 


For 0 = 0 it is reducible and can be reduced to the implicit Euler method 

with r(0) = -1. For 0=1 it is reducible too and the reduced method is the 

explicit Euler method with r(l)» 1. For 0€ (0,1) one finds r(0)= 1/0. 

In particular, for the trapezoidal rule, where 0=1/2, it follows that 

r(l/2) = 2. This result is in agreement with the fact that the trapezoidal 

rule is not B-stable, see [12]. To compute the stability region we observe 

that K(yl) = (l +y0)/(l - (l-0)p). Hence S = D(r (0)) with r (0)= 1/(20-l). 

s s 

Therefore one has 


r(o) «-1 «r (0) 

8 

D(r(0)) is a proper subset 
of D(rg(0)) 

r(l ) = 1 « r^(l ) 


implicit Euler, 
for 0 < 0 < 1 
explicit Euler . 


Note, however, that if we define, 

y - y - (l-0)hf(t ,y ), 

then satisfies the "one-leg” difference equation, 

Cl • 

see [U], It is well known [U] that this one-leg method is A-contractive 

(B-stable) when 0< 6< If 2 is defined a?:alogously to 5^ , it follows 

£ n n 

that 

"Cl - Cl" ^ 

The conclusion is that, for 0<6<|, the e-method, although it is not B-stable, 
is A-contractive with respect to a different problon-dependent , metric. 

Example 4 . The most general two stage second order explicit Runge-Kutta 
method is characterized by 

^“(l/ 2 o o)* b’’-(l-a.«). «*0, 


see [6, p. 121]. If a* 1 the method is reducible and thus not circle 
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contractive. However* for a€ (0*1) one finds by an easy calculation that 

(5.1) r(a) ■ 2/(1 ] -, 7-3) . a€(0*l). 

Here r(o) depends truly on a* and r(a)< 1 for ^is is in contrast 

to the stability region S which is independent of a. In fact 
S»{y€C| |l + y + y^| <1} and thus 

(5.2) r (a) ■ 1 for all o*0. 

8 

It is well known that S is bounded for explicit methods. Hence r is i>ositive 
for explicit circle contractive methods. How large can r actually be? 


Theorem 5.1 . Assume an explicit m-stage Rimge-Kutta method is r-circle 
contractive. Then 

(5.3) r < m . 


Moreover* equality is only attained if 


(5.U) K(ul) - (1+J)“. 

xn 

which implies that the error order is one. The method with 

bj^ ■ 1/m i ■ 1 *2*. . . *m. 


(5.5) 



for i < j 
for i > j 


attains equality in (5.3). 


Proof. In [8] it is shown that r_sm with r_*m if and only if (5.*») holds. 
Thus from r<r^ follows (5.3) and (5.^). If a Runge-Kutta method has error 
order p, then K(wl)-e'^ ■ 0(y^^^). For the specieil K(yJl) of (5.**) we find 
e*^-K(yJl) ■ ^^(w®) and thus by (2.U) p»1. An easy calculation 

shows that B ^QB ^ ■ -^I for the method given by (5.5). Thus by (3.26) 

SI SI 

one has r«m and equality in (5.3) holds. 


Let us now consider the same problem for implicit methods. Burrage and 
Butcher [1] have investigated algebraic stability and shown that there are 
implicit m-stage Runge-Kutta methods of order 2m* 2m- 1 and 2m-2 which are 
algebraically stable* that is r is nonpositive. The following theorem 
gives a relation between the size of a negative r and the accuracy of tha 
method. 


Theorero S.2 . Assume the Ruage~Kutta method is r-circle contractive with 
r < 0 and 


(5.6) 

Then 




p ■ 1 , c < 0 
and 

r < 1/2c . 


Proof . Let R be the radius of curvature of dS at u^O. Since D(r)cS one 
has that 0 < R S -r. It remains to be shown that 


(5.7) 


if p> 1 


if P-1. 


Lot dS be given in a neighborhood of 0 by the equation w * c(t) + it. 

C(t.) is implicitly defined by |K((c(t) + it)jQ| ^ ■ 1. Using (5.6) we find 


(5.8) 


1 . 

(l ♦ c(s(t ) - It)^’ ♦ 0{« (t ) - it)'*®) 

Implicit differentiation of (5.8) gives 

V(0) - 0 

C«(0) - / ° if P> 1 
l-2cifp«1 

and hence (5.7) follows ismediately. 


Note that for implicit r-circle contractive methods with a nonpositive r 
the absolute value of r increases as the accuracy increues. 
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6. Calculation of r for some explicit methods 


We omit the algehredcally stable methods given in [1] and restrict ourselves 
to the explicit methods listed in [9]« 

All aeoond order two etage methods are contained in Example 4. 

Third order formulae. Note that for all these formulas, rg*1.25. 

Classic form 


A 


b 


T 





r«0.5 


Nystrom form 


Heun form 



- ( 'C ^ rwO.927 



iK 0 % ) 


This method has b 2 ■ 0 and is irreducible. Thus it is not circle contractive. 


Reilston's optimum third-order form 

/O 0 0 \ 

A - I 0 0 ) 

\0 X 0 / 

b^-^(2 3 4 ) r« 0.099 

Kuntzmann's optimum third-order form 

/ 0 0 0 

A>( 0.4648162 0 0 

\-0. 0581 020 0.8256939 0 

b’^ - ( 0.2071768 0.3585646 0.4342585) r«0.847 

Fourth order formulae. Note that for all these formulas, r ~1.4. 

B 


/O 0 0 0 \ 
^ 0 0 0 ^ 
V 0 0 0 y 

^ 0010 ' 


r« 1 


Classical form 




' 20 ' 


Kutta fora 


Gill form 


/O 0 0 Ov 

(4 ? S 0 

' 1-1 1 


^ ^ ^ H ^ ) ra. 0 .U 61 . 


0 0 

^ 0 

,(\/J-l)/2 (2-\/5)/2 

0 -\/?/2 
^ (2-V5)/6 (2+V5)/6 



Kuntzman optimum fo\irth order form 



- 0 

0 

0 

J_ / 

88 

0 

0 

220 V 

-33 

165 

0 

\ 

95 

-75 

200 

1 / 
3^ ^ 

55, 

125. 

125 


) 


r«0.586 


re* 0.698 


Ralston's optimum fourth order form given in [9» p.38] is not circle 
contractive since b2~ -0.55198066 < 0. 

Concerning methods of order exceeding four* see Theorem U.U. 

Finally, we recall the remark, made in Example 3 of Section 5> One can 
perhaps find a larger value of r with a different metric. Therefore, our 
values of r must not be considered as a final verdict in the comparison 
of methods. Our conditions are sufficient for good behaviour on certain 
non-linear problems, rather than necessary. 


It is also possible that the picture can be brighter for some methods, if 
we relax our requirements a little in other respects, e.g. by practically 
reasonable regularity assusq;>tions for the function f. One can perhaps 
"break the barrier" expressed in Theorem by small changes of our 
definitions. 


Hyman [7l has recently reported some interesting empirical evidence of the 
shortcoeiinfs of the linear stability theory as a guide-line for the 
behavioisr of Runge-Kutta methods on non-linear problems. We have not yet 
had the opportunity to study his results from our point of view. 

More theoretical and experimental research is therefore needed to test the 
practical relevance of our analysis. 
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